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1.  INTRODUCTION 

The  system  for  quantum  computation(QC)  that  we  are  proposing  is  modeled  on  the 
ion  trap  system  for  QC1  where  the  ion  qubits  are  coupled  among  themselves  through  their 
collective  oscillations  in  the  ion  trap,  which  plays  the  role  of  a  resonant  bus.  The  bus  for 
the  proposed  Josephson  junction(JJ)  prism  qubit  QC  is  a  single  mode  LC  resonant 
loop2'3'4'  Magnetic  flux  from  the  loop  threading  a  qubit  will  produce  a  loop-qubit 
interaction  via  the  Bohm  Aharonov  effect.  At  their  preferred  operating  point,  the 
proposed  qubits  have  negligible  circulating  persistent  currents  and  are  kept  by  symmetry 
from  interacting  with  their  coupled  measurement  SQUIDs.  By  pairing  the  physical 
qubits  into  logical  qubits,  it  is  possible  to  maintain  the  logical  qubits  in  a  decoherence 
free  subspace(DFS)5,6  7  that  nulls  the  decoherence  effects  of  uniform  external  flux 
perturbations,  and  the  perturbations  associated  with  computational  Molmer 
Sorensen(MS)  bichromatic  gates8'910.  The  only  time  a  logical  qubit  should  interact  with 
the  environment  is  during  the  unitary  gate  rotation  that  initialize  the  qubits,  and  during 
the  final  Hadamard  rotation,  which  allows  the  qubit  to  interact  with  the  measurement 
SQUID  gradiometer. 


2.  JOSEPHSON  JUNCTION  TRIANGULAR  PRISM  QUBITS 


The  physical  qubits  that  we  shall  consider  are  based  on  a  simplified  version  of 
qubits  considered  in  a  previous  paper4.  We  assume  that  the  junction  critical  currents  and 
the  circuit  dimensions  are  small  enough  so  that  the  dimensionless  self  and  mutual 
inductances  of  the  qubit  circuit  in  units  of  Lc  =  <D0 l{2nlc)  are  small  and  can  be 
neglected11'12.  The  Hamiltonian  for  the  qubit  shown  in  Fig.  1  may  be  written  in  term  of 
the  junction  gauge  invariant  phases  Cl ,  critical  currents  J a ,  and  capacitances  C,  as 

//  =  X[y<£2(®0/2.7r)2/Ci -EJ(Ja/Jc)cos(^j)^  where  the  flux  quantum  and  junction 

characteristic  energy  are  written  as  <D0  =h/2e  and  E ,  =  ./  / 2n  . 

The  phase  constraint  equations  for  the  two  cells  are  given  by  <£j  -  £2  -  =  <fa x  and 

£2  -  £j  +  =  fa2  ■  Defining  the  symmetric  and  antisymmetric  combinations  of  the 

external  magnetic  fluxes  and  phases  as  &  =  (£l+£1)/2,  y/  =  (£,  - £,)/2  and 

fa  =  Wax  +  fai\  fa  =  Wax  ~  fai)  >  with  <t,  =  2.n($?i /O  0)  yields  ^=2y/- </>a  -  <j>s  and 
^4=2 y/-<ta+<ts..  We  allow  for  the  possibility  of  differing  junction  areas  by  writing 
JCi2=J±AJ  ,JC3A  =  j  ±  Aj ,  Cj  2  =  C±  AC ,  and  C34  =C±AC. 

In  order  to  obtain  the  unperturbed  Hamiltonian  and  its  wave  function,  we  assume 
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Figure  1.  a)  JJ  prism  qubit  with  junction  critical  currents  Ja,  gauge  invariant  phases  £  and  external  magnetic 
fluxes  <l>i  indicated,  b)  Resonant  LC  bus  symmetrically  coupled  to  qubit  pairs,  c)  Antisymmetric  coupling  of 
adjacent  buses  to  a  logical  transfer  qubit. 


equal  area  junctions  and  circuit  loops,  zero  external  antisymmetric  flux,  and  a  constant 
symmetric  flux  given  by  <D°  =  (</>°  /2u)<t>0  .  The  Hamiltonian  may  be  written  using 

9  and  y/  variables  as 


H0{9,y/)  =  \Med2  +\Mvy/2  -  Ej [2 cos($) cos( y/)  +  2(j  /  j)cos(</>°)cos(2ys  -  <j>a )]  (1) 

where  we  have  defined  M¥  =  2C(l  +  4/?)(cD0 Hn^)  =  (l  +  4p)tr/4Ec ,  p  =  C/C, 

Me  =  2 C  (<D0 /2;zj‘  =fi2/  4 Ec ,  with  the  capacitive  energy  Ec  given  by  Ec  =  e2  /2C  .  A 
contour  plot  of  the  potential  energy  V(0,y/)  =  -2cos(0)cos(y/)  +  r  cos(2y/ -  </>ci)  is  shown 
in  Fig.  (2a)  for  r-  0.675  where  r  =  -2(j /  J)cos(^s°) .  Approximate  eigensolutions  are 
obtained  by  making  a  Hartree  approximation,  where  the  full  two  dimensional  wave 
function  is  approximated  by  a  separable  wave  function  KY(d.t//)  ~  Q(d)kY(t//) .  Assuming 
that  the  motion  in  the  9  direction  is  in  its  ground  state,  the  cos($)  term  in  V{9,y/)  may 
be  replaced  by  its  expectation  value  r/H  =  (Qe  |cos(^)|0#) ,  leading  to  an  effective  qubit 
potential  in  y/  given  by  V0{y/)  =  -2r]Hcos{y/)  +  >~cos{2y/) .  Writing  9  and  y/  in  terms  of 
the  momentum  operators  M g9  =  pe  =-ihd / 89  ,  M  y/  =  pw  =  -ihd/dys  with 
m¥  =  M v Ej  /  tr  ,  the  unperturbed  effective  ID  Hamiltonian  leads  to  the  time  independent 
Schrodinger  equation 

-d2xV(y/)l dy/2  -2m  [2rjH  co^{y/)-rcos{2y/)^¥{y/)  =  2my/{Eil EJ)xY{y/)  (2) 

which  is  in  the  form  of  the  Whittaker  Hill  equation.  Solutions  for  the  eigen  functions  and 
energies  for  the  Whittaker  Hill  equation  have  been  found  by  Unwin  and  Arscott  terms  of 
continued  fractions13. 
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Figure  2.  a)  Contour  plot  of  V(0,y/)  for  r  =  .675  b)  Energy  levels  for  V(0,y/)  in  units  of  Ej 


3.  COUPLING  THE  RESONANT  LC  BUS  TO  THE  QUBITS 

The  qubits  are  coupled  among  themselves  via  their  (symmetric)  mutual  inductance 
coupling  with  a  resonant  LC  bus  as  shown  in  Fig.  (lb).  The  Hamiltonian  for  a  bus 
oscillator  with  inductance  L  and  capacitance  C  may  be  written  in  terms  of  the  phase 
operators  cps  and  p,.  as  i/“c/(<t>0/2^)2  =p]l2M(p  +  M va\c(p2 1 2  with  =  C(<1>0  Htt)2  , 

coLC  =1/  ~JLC  ,  and  M  q>s  =  ps  =  —id/d<ps .  In  terms  of  annihilation  and  creation  operators 
<ps  =  ^Til  iMp&jf.  ( a  +  at)  and  ps  =  i^Mj\(DLC  I2(a1  -  a) .  In  order  to  carry  out  unitary 
qubit  gates  on  the  ith  qubit,  we  include  an  external  microwave  pulse  8<j)^\t)  to  be  added 
to  as  well  as  an  external  microwave  or  DC  pulse  S(f>[‘lxt(t)  to  be  added  to  (j>l')0 

along  with  ku)(ps ,  the  portion  of  the  bus  flux  that  interacts  with  the  ith  qubit,  where 
£(,)  =  (Z(,) /Z)  with  li"  the  inductance  of  the  ith  coupling  loop.  For  the  case  with  no 
defects,  we  may  write  the  Hamiltonian  for  the  ith  perturbed  qubit,  with  the  definitions 
coo  =(0|cos2^|0),  cn  =  (l|cos2^|l),  501  =  (0|sin2^|l) ,  £\  =  {j\H0  \j)l Ej  as 


(0 


1 bus  0  1  (,)(ZZ'°  +  H™  +  H osc) / E jl(n  0  lhus 

S  |(e0  +  eji'0  0  lbus  +  -  ex )df  ’  0  lbus  +  1  bushcoLC  IEJ(a]aj  +  |)ltal  v,  x 

+ j  pMg  i(,)  0 1'”“  + \  PMg  \bus  s$l)2  ii!,s  0  i(0 

+ih"{iKc„  0  +  cn)[(/0(^°)  -  1)  +  1)  -  tan  (f>l  J,(Sf)  sin^f]}^  0  l(i) 


+  1A"S  { j  7-(c00  -  cn )[  ( J0  (SI1  > )  - 1)  +  (J0  (sn  -  1)  -  tan  fs  J0  (%'>)  sin  S^>  ]  }1 


-(On  litas 


1(0 


+{p£°  /(l  +  4p)2EJ}(^0(Vfl)\d/dVft 

+  iA">s01(cos S</>^  sin 5(j>(p  -  tan </>°  sin  SfJ]  sin<%(0)}i6'“  0  <7<° 


(3) 
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3.1  Logical  Qubits  for  a  Decoherence  Free  Subspace 

The  Rabi  frequency  Q1'1  for  coupling  a  bus  excitation  to  the  ith  qubit  via  the  interaction 
term  smS^’  sincf^01  a^1  — »  S^1’  cp^  '  is  given  by 

Q(,)  =  [rtan(^° )(?]'’  s10^JTuolcL/2  I2L)] .  Since  Q(,)  oc  1/ 4n  ,  to  increase  Q(,)  fora 
particular  value  of  N  it  is  necessary  to  increase  either  S(f>l‘)  or  Z(,)  or  both.  Increasing 
L(,,(by  forming  multilevel  spiral  inductors)  is  only  mildly  effective  due  to  the  (Zl,,)T 
dependence,  while  increasing  the  amplitude  of  <%(l> ,  while  effective  has  the  drawback 
that  the  Bessel  expansion  term  r^O|cre2^|0^l|cos2^|l^((/0(<5(',)-l)(5{')  =  E^j  of1 

is  excited  concomitantly,  and  will  lead  to  phase  errors  when  S1^  is  not  small. 

This  type  of  phase  error  (as  well  as  any  symmetric  external  flux  fluctuations)  can  be 
nulled  by  using  encoded  logical  product  qubits  with 
{|0z)  =  |0l)  =  |0)<8>|l),  |li)  =  |l0)  =  |l)®|0)}  or  by  employing  the  two  entangled  Bell 

states  {|6i)  =  (|0l)  +  |l0))/V2,  |l^  =  (|0l)-|l0))/V2}  as  logical  qubits.  For  either 

type,  applying  and  to  the  first  (i)  and  second  (i+1)  physical  qubits  of  a 

logical  qubit,  will  result  in  a  positive  phase  shift  in  one  physical  qubit  canceling  a 
negative  phase  shift  in  the  other  if  their  two  Rabi  frequencies  (XT1'1  and  (XTl‘+1)  are  equal. 

It  can  be  shown  that  this  type  of  cancellation  will  occur  during  MS  gates  as  well. 

The  Bell  state  qubits  have  the  property  that 

K  =(5,  +^2)|li)  =($  -£-„)(I) -te-Sof  and  (0L  +^2,|0i)  =  (l,  |  Hm  +//f2)|li)  , 

i.e.  their  eigen  energies  are  exactly  equal  (if  they  are  properly  formed).  However  if  the 
junction  or  flux  properties  of  the  two  physical  qubits  differ,  the  transition  matrix  element 
tr  will  be  non  zero  and  there  will  be  tunneling  between  the  two  states.  Transforming  to  a 
basis  in  which  tr  =  0 ,  yields  the  product  states  with  zero  tunneling  but 
(0L  |//(1)  +//,2)  |0Z)- (lz  |//(1)  +H{2)  |lL)  =  (sl -s0)(2) -(£,  -  ff,,)11’  so  there  may  be  non 
zero  time  evolution. 

3.2  Molmer-  Sorensen  Gate 

Sorensen  and  Molmer  have  given  a  detailed  description  of  a  bichromatic  excitation 
scheme  that  enables  multiple  qubit  gate  operations  to  be  carried  out  for  ion  trap  quantum 
computers8’910.  Their  scheme  is  desirable  since  it  allows  for  relatively  rapid  gates  to  be 
carried  out  via  an  oscillator  bus,  while  returning  the  bus  to  its  pre-gate  initial  configuration 
irrespective  of  the  initial  state  of  the  bus.  For  the  case  of  logical  JJ  prism  qubits  coupled  by  a 
resonant  LC  bus,  the  MS  scheme  is  desirable  since  it  allows  logical  qubit  gates  to  be  carried 
out  without  exciting  figure  8  currents.  As  the  gate  operators  for  the  ion-trap  case  have  been 
set  forth7”'  we  display  the  needed  gate  operators  for  the  slightly  different  case  of  the 

interaction  Hamiltonian  sin(5$,)sin(5$',a1r'*  while  relaxing  the  condition  that  all  qubit  interactions 
have  the  same  Rabi  frequency.  Expanding  the  sin<5i^,)  and  sin  St/>^  terms  to  lowest  order  in 
St fP  and  8$*  =4 / L)  and  making  the  rotating  wave  approximation  leads  to  an 
interaction  term  V(t)  =  -^2Q/t(^)[^(cos(®iC -5)t  +  ps(sin(®//;  -  J)t]  =  Jx(jj>)[f{t)(ps  +g{t)ps\ 


5 


where  Jx($)  =  \'n.d*&+  +  e  ,<4<7i](Q,./Q)  withfl  =  [rtan(f^°)<^0 sl0^tia>LCL/ 2  ju(Lsile / 21)] 

i 

and  Q  being  the  individual  and  average  Rabi  frequencies  of  the  interacting  qubits.  For 
the  MS  bichromatic  excitation,  is  made  up  of  signals  from  two  microwave  sources 
with  source  phases  </>0)  and  frequencies  co±  =  cul0±S,  where  co10  =  (s,  -  £0) /  ft  and 
<?□  coLC  ±20..  In  Ref.  10,  the  exact  evolution  operator  for  V(t)  is  shown  to  be  given 
exactly  by 


U (t)  =  exp (~iA(t)J2x )  exp (~iF(t)Jx  x)exp (~iG(t)Jx  p )  (4) 

with  F(t)  =  \'0f(T)dT  =  -'j2Osm[(coLC -S)t\/(wLC -5),  G(t)  —  \'0g{r)dr  and 
A(t )  =  -\'nF(r)g(T)dT.  For  a  pulse  length  given  by  ( a>LC  -S)tK  =  2nK  where  K  is  an 
integer,  F(tK )  =  G(tK )  =  0  and  the  bus  can  be  returned  to  its  initial  vibrational  state. 
Choosing  values  for  K,  8,  and  FI ,  yields  A(tK)  =  -0227rK /(coLC  - S)2  and  conversely 
tK  =  yj- A2xK  / O  .  The  single  (logical)  qubit  gates  can  be  composed  as 

X  =  2  J2M>  =  0)=  HA2  +  a±)\  + 

Y =  2Jl (*  =  0,  <4  =  \n)  =  }(A?  +  A2)Y -  A 

where  At  =  Q  /  O  .  Operators  for  product  and  Bell  logical  qubits  are  given  respectively 
by  {o',,  =cr'1'o-‘2),  <t;>.  =ct'V'2),  alz=-i(7lx(7r  y  }  and 

{°Lx  =  -i dL,  °Ly  =  ~<*(y  ’>  ^Lz  =  ¥  or  the  product  logical  qubits,  the 

CZ  gate  may  be  written  as  U^U)  =  V7 exp(icr^)cr^)^-/4)exp(-;cr^V/ 4)exp(-/cr^)^/4) 
where  the  term  exp(/cr|/r)<T^fV/4)  can  be  expressed  in  terms  of  a  four  (physical)  qubit 
Molmer-  Sorensen  gate  XX  =  2  J2x(j>x  2  3  4  =  0) 

Ucz'n  =  exp(i>/ 4)exp(iXf^-/4)exp(;cr^V/ 4)exp(/cr}“V/ 4)  (6) 

When  the  Rabi  frequencies  for  the  four  physical  qubits  are  not 
identical,  mixed  terms  like  cr[a4  will  arise  with  non  zero  coefficients  proportional  to 
yl  =  JLr  -  I .  If  the  Rabi  frequencies  for  each  physical  qubit  of  a  logical  qubit  pair  can  be 
made  equal,  then  FI  can  be  defined  as  the  average  Rabi  frequency  with 
Y  =  Yi  =  Y2-  ~7i  =  ~Ya  ■  Terms  that  are  proportional  to  Yi+Y2+ 7 3  +  7a  w*b  be  nulled, 
and  the  remaining  terms  will  be  proportional  to  either  a[l)a[2) ,  ax"a[4> ,  or 
cr^’cr^’cr^cr*41  that  can  be  accommodated  by  small  changes  in  pulse  lengths.  If  the  two 
Rabi  frequencies  (XT1' *  and  <XTt,+1)  for  each  logical  qubit  pair  are  equal,  as  is  required  to 
cancel  out  SOU)o’[')  and  $Q<'+1,c7<:'+1  phase  errors  generated  during  a  gate,  then 
yfcos  (<0(4>  -  cf!,)[<OT  =  j™  cos(<0(c<2)  -  c^ns^f  is  satisfied.  If 
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y'f)sin(^J(g))Sg])[<^1)]=  y'2)  sin(^))5jf  [<?|2)]  is  satisfied,  then  the  Rabi  frequencies  for 
qubits  1  and  2  will  be  equal.  For  both  constraints  to  hold  concurrently,  the  equality  of 
Eq.  (7)  must  be  satisfied,  which  is  possible  to  achieve  by  varying  4  (or  SR )  if  (f)('(j  and 


(4  IS?)  =  tan(^>)^(C  -  C)/ tan(^)<(42)  -  c[?)  =  Rr <V(2)) 


(7) 


(Zi^’are  fixed  ,  as  shown  in  Fig.  (3). 
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Figure  3.  a)  y  Transition  frequency  (for  E j  /  h  =  200GHz  )  and  const.xRabi  frequency  Q  vs.  r. 
b)  as  a  function  of  for  r(I)=0.60,  0.65,  and  0.75. 


4.  CIRCULATING  CURRENT  PATTERNS  FOR  THE  QUBIT  STATES 

The  circulating  current  patterns  associated  with  the  |0^  and  |l)  states  can  be 
derived  by  using  <^3  4  =  ly/  -  (j)a  +  tj>s  and  by  assuming  <j> to  be  a  small  perturbation, 
which  then  yields 

4  =(°  |  Ria  sin  4  3,4 10) 

4  =  (ikca^in^ll) 

with  k0  =  sin^°  cosfzi0  (0|cos2^|0^  and  kx  =  sin^°  cos^l7l|cos2(//|l^  .  In  the  Flartree 
approximation  if  AJ  =  A/  =  0  the  expectation  value  of  JX=JCX  sin  4'  in  the  |l^  or|0^ 
state  of  y/  is  equal  to  zero.  The  circulating  current  patterns  associated  with  the 
unperturbed  |0)  and  |l)  states  are  thus  clockwise  currents  of  magnitude 

4  =  -4  =  ~jK 0  and  4  -  -4  -  ~JK\  around  the  perimeter  of  the  qubit,  and  are 
plotted  vs.  r  in  Fig.  7b. 

Since  the  probability  density  for  the  two  superposition  states  |±)  =  4^(|0)±  |l))  is 

centered  around  the  potential  well  minima  at  y/  =  +(//  ,  circulating  currents  for  these 
states  are  close  to  the  classical  equilibrium  currents.  Writing 

Km  =(0|sin2(//|l)cos((4°  -$),  circulating  currents  for  the  |±)  states  may  be  reexpressed  as 


=  +  y'AT0 ( 1  ±  cos <fP  sin (j)a  / sin^s°  cos^u ) 
=  +  /at,  (1  ±  cos  $  sin  (j)a  /  sin  fzi”  cos  <j>a ) 


(8) 
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4+>  =  (+|^C3  si<3  |+)  =  ~\j(K 0  +  KX  )  +j %> 

4+>  =  (+1^4  sin  C|  ■ +)  =  +i  M)  +  *q )  +  j. Km, 

The  circulating  current  is  thus  comprised  of  a  perimeter  current  |y(/^+/q),  which  has  the 
same  sign  for  the  |d^  states,  and  a  figure  eight  current  jxm  that  changes  sign  for  |±^  states. 

4.1  Initializing  the  Qubits  and  the  Effect  of  Critical  Current  and  Geometric  Defects 

All  of  the  qubits  are  initially  allowed  to  relax  to  their  respective  ground  states.  For 
a  system  in  which  all  of  the  cell  and  junction  parameters  are  identical,  every  even 
numbered  qubit  is  given  an  x  rotation  of  180°.  In  order  to  avoid  phase  errors  during  these 
initializing  single  qubit  gates,  the  gates  should  be  carried  out  with  long,  small  amplitude 
microwave  pulses.  Each  {iodd,  iodd+i}  pair  of  qubits  is  then  designated  as  a  product  logical 
qubit  {I, II, III..}.  To  initialize  the  Bell  state  logical  qubits  a  further  logical  (MS) 
Hadamard  gate  is  needed  for  each  logical  qubit.  A  test  measurement  should  show  that 
there  is  no  antisymmetric  flux  coupled  to  the  SQUID  gradiometer  detectors 

For  a  qubit  in  which  there  is  some  combination  of  junction  critical  current  defects 
and  cell  geometry  defects,  a  term  in  the  Hamiltonian  will  be  generated  that  is 
antisymmetric  in  y/  and  6 

AH  =  -{2/sin </>rlgeam  cos $  +  2A/cos^  ,,fM,„  sin^°}sin2(//  -  2 A/ sin <9 sin;//  ^ 

=  -Kdefect  sin  2y/  -  2  A  J  sin  9  sin  y/  =  A  V“  +  A  Vb 


4_>  =  ("I  JC3 sin  £  | ' ~)  = ' ~\j(Ko  +  *i)  -  J  Km 
4_>  =  (-kc4  si<4|-)  =  +\AK0+K\)~jK0l 


The  effect  of  AVa  on  the  potential  is  to  tilt  it  in  the  (//direction  as  shown  in  Fig.  4b, 

while  the  effect  of  AF*  is  to  skew  the  potential  as  shown  in  Fig.  4c.  For  AVh  an 
accurate  numerical  solution  for  the  2D  potential  has  been  undertaken  but  not  yet 
completed,  we  will  thus  approximate  the  perturbed  wavetlmctions  using  first  order 
perturbation  theory. 


Figure  4a.  a)  Symmetric}  kq  ,  }  and  antisymmetric  { /u01  }  current  amplitudes  vs.  r. 

b)  Contour  plot  of  V(0,y/)  for  AVa  ^  0  .  c)  Contour  plot  of  V{0,i//)  for  AVh  ^  0  . 


In  a  Hartree  approximation  there  are  four  types  of  states  when  labeled  by  parity  with 
respect  to  9  and  y/  ;  \m,n)=  {\even,evenj ,  |  even,odd^ ,  |  odd, even),  and  |  odd, odd}} 
where  m  and  n  are  the  number  of  nodes  in  9  and  (//  .  The  only  type  of  (separable)  states 
contributing  to  the  summation  that  will  survive  in  a  first  order  perturbation 


expansion  j  n)  =  |m°\+  X  (*°  I A  Vb  |n°)/(£°  ~  E°k)  for  |o)  and  |l),  are  |  odd, odd)  and 

1  '  k*n  '  1  ' 

| odd, even')  type  states  respectively,  where  |/70^  are  unperturbed  states.  Keeping  only 
the  lowest  energy  terms  1 2° ^  =  1 1, 1^  and  1 3°^  =  1 1,0^ ,  an  approximate  solution  may  be  written 
as|o)  =  |o°^  +  a|2°^and  |l)  =  |l0^  +  /?|30^  where  a  =  2AJ  (2°  |sin0sin(//jo°y(£jo) -E^)  and 
P  =  2A J  ^3°  | sin 6> sin ^ 1 1° ^ /(i?30)  -  E ^0>) .  The  matrix  elements  of  AVh  are  then  given  by 
(0|AK* |0)  =  2«(0° | AK* |2°),  (l|AFi’|l)  =  2^(l0|AF'’|3°),  and  (0|AF*|l)  =  0.  Critical 

current  defects  injunctions  1  and  2  will  thus  lead  to  a  second  order  (in  A J / J  )  correction 
in  the  energies  of  the  two  qubit  states,  but  it  will  not  lead  to  transitions  between  them. 

The  most  direct  way  of  nulling  the  A  V"  term  is  to  add  a  small  external  antisymmetric 
‘trim’  flux  to  the  existing  flux  so  that  the  resultant  antisymmetric  flux  <f =<j>agmm  +  </>a  lnm 

will  satisfy  tan^°  =  -(A // /)  tan  (f,  .  The  net  result  of  having  (zSj01  determined  by  this 
equation,  is  that  small  figure  eight  currents  that  were  present  initially  in  the  |0^  and  |l) 
states  and  could  interact  with  the  SQUID  detector  are  now  nulled,  while  perimeter 
currents  in  the  |0^  and  |l^  states  will  be  slightly  decreased  by  the  last  factor  in  Eq.  (8) 

which  simplifies  to  (1  ±  A/7y)(l  +  cos^s°  sin^a /sin^“  cos^a)  =  [1  -  (A//y)2] .  If  A JIJ  is 
not  zero,  then  there  will  be  a  multiplicative  factor  of  (1  +  cTc,,  /  c00)  for  the  |0^  state  and 
(1  +  /?2c00/cu)  for  the  |l)  state. 

If  the  trim  current  is  not  applied,  then  there  will  be  a  non  zero  figure  eight  current  of 
magnitude  kq  ( Aj + j  tan  </>a  gixm  /  tan  (j)] )  for  the  |0)  state  and  K]  ( Aj  +  j  tan  </>a  geom  /  tan  </>’ )  for  the 

l)  state.  From  Fig.  4a  it  can  be  seen  that  if  the  qubit  is  operated  around  r*,  then  k0  and 
kx  are  small  (~  equal  and  opposite)  and  0.Rabi  is  near  its  maximum  at  r  ~  0.75  . 

The  Flamiltonian  H0  +  AV  can  be  diagonalized  by  a  unitary  rotation  to  a  new  basis. 
Since  the  unperturbed  wave  functions  for  the  two  basic  states  x¥0(G,y/)  =  {6,y/\0^  and 
=  (0,i// |l)  are  both  symmetric  in  9 ,  the  matrix  element  of  the  last  term  in  AV 
will  be  zero.  Only  junction  defects  Aj  ,  associated  with  outer  junctions  3  and  4,  and  a 
small  antisymmetric  flux  (j)a  due  to  different  cell  sizes,  will  then  affect  the  operation  of 
the  qubit. 

4.2  Connecting  Buses  into  a  Network 

It  is  possible  to  couple  groups  of  buses  in  an  open  branching  network  using  transfer 
qubits  as  first  described  in  Ref.  (4).  A  diagram  of  one  possible  coupling  scheme  is  shown 
in  Fig.  lc  where  bus  1  is  coupled  symmetrically  to  the  transfer  qubit  and  bus  2  is  coupled 
antisymmetrically  in  order  to  prevent  direct  mutual  inductance  coupling  from  bus  1  to 
bus  2.  Using  U^ju>  from  Eq.  (6),  an  expression  for  U\jjk>  may  be  written  using  the 

Fladamard  operator  Hk  as  U(jNk)  =  H kU{j7k) H k .  Writing  Hk  =  -i(f‘ul4)<T'y  e'u  2)<7'f  ,  the 

CN  operator  may  be  expressed  as  U(cJNk)  =  e  a,r,4)CTlr  which  requires  4 

(bichromatic)  microwave  pulses.  Since  {x(l),  p(l)}  commute  with  {x(2),  p(2)} ,  the 
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expression  for  Aft)  in  Eq.  4.  is  replaced  by  A(t) +  A(2)(t) .  For  both  bus 
oscillators  to  return  to  their  pre  gate  values  both  of  the  constraints 
(co(X)lc— S)tK  =  Km2jr  and  {co^LC  —  S)tK  =  K(1) 2/r  with  integer  K (/)  must  be 

satisfied  simultaneously.  For  K(2)  —  -K{1)  a  solution  is  S  -  (a>(i)LC  +  a>a)LC) /  2 , 
tK  =  AkK(X)  /(o)(l)LC -ty(2)ic) ,  which  yields  An}(t)  = —A{l)(t)(-l)P  where  P  is  the 
parity  of  the  bus  l-bus2  coupling.  Thus  U(tK)  =  exp(-2iA(lAtK)a{xl)(T(x2)) .  . 

For  the  2  logical  qubit  CZ  gate,  if  one  of  the  logical  qubits  is  a  transfer  qubit,  its 
coupling  to  a  second  bus  will  engender  a  single  logical  qubit  rotation  that  is  in  the  proper 
direction  for  the  CZ  gate.  Thus  if  logical  qubit  3-4  is  a  transfer  qubit  connected  to  a 
second  bus,  the  expression exp(i7tXX / 4)  =  exp(-in !  A)exp{-i7t<j{))<j^)<j^)<j\*)  / 4)  picks 
up  a  factor  exp[i7ra<x3)aix4>  /  4]  obviating  the  necessity  a  separate  gate.  If  there  are  two 
transfer  qubits  involved,  no  single  logical  qubit  gates  are  needed  at  the  transfer  qubit 
sites.  This  reduces  the  number  of  separate  microwave  pulses  for  a  CN  gate  with  one  bus 
transfer  from  5  to  3  and  also  reduces  the  time  needed  to  carry  out  the  gate  by  2 tK  . 

To  carry  out  a  CN  gate  over  a  number  of  buses,  it  is  necessary  to  first  carry  out  a 
sequence  of  Uc'Nn ’s  starting  from  the  control  qubit,  through  all  of  the  intervening  transfer 
qubits  to  the  target  qubit,  followed  by  a  final  sequence  of  U(ff, )  from  the  original  qubit 
to  the  last  transfer  qubit  in  order  to  reset  the  transfer  qubits  to  their  original  ground  states. 
The  sequence  involving  2  transfer  qubits  would  thus  be  1  between  (  /',?,),( tI,t2),{  t2,k) 
followed  by  resetting  transfer  qubits  with  1  applied  to  (j,tf)  and  (th  t2).  These 
maneuvers  may  be  chained  across  as  many  coupled  buses  as  required.  The  wave 
equation  for  the  entire  system  may  be  written  as  =  X  „c„  j'P" )  with 

| 'P”  )  =  nil  |  Xa’ )  ®  n  -:,  |  X^ ^  ®n&  I X(k\n  )  >  Where  the  qubits  and  bus  are  restricted  to 
the  two  lowest  energy  states  al  =  {0, 1 } ,  =  {0, 1 } ,  yk  =  {0, 1 }  and  where 

«={«i  a2  ....aN  P2,...pN  ,yx,y2,....ymJ  with  |  j0)  =  |6Z),  |i,)  =  |n),  |j0)  =  |0flUI)» 

and  |  X\ )  =  |  )  •  Because  the  MS  gates  return  the  buses  and  (along  with  resetting)  the 

transfer  qubits  to  the  initial  state  they  were  in  before  the  U^k)  gate,  and  the  logical 
qubits  do  not  evolve  in  time,  it  should  be  possible  to  maintain  quantum  coherence  across 
the  set  of  coupled  buses.  This  allows  open  branching  networks  of  coupled  buses  to  be 
formed,  enlarging  the  possible  number  of  interacting  qubits  in  a  single  quantum 
computer. 

CONCLUSIONS: 

We  have  shown  that  it  is  possible  to  design  a  quantum  computer  comprised  of 
Josephson  junction  prism  qubits  that  are  coupled  via  the  flux  of  a  resonant  LC  bus.  that  is 
capable  of  carrying  out  quantum  computations  in  a  decoherence  free  subspace.  The 
effective  Flamiltonian  is  zero  between  gates,  so  there  is  no  dynamical  evolution  of  the 
system.  The  |0)  and  |1)  physical  qubit  states  are  characterized  by  circulating  currents  that 
follow  the  perimeter  of  a  two  cell  qubit,  and  are  thus  symmetric  with  respect  to  the  flux 
threading  the  cells  of  the  qubit.  The  logical  qubit  states  are  made  up  of  linear 
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combinations  of  |0l)  =  |0)<8>|l)  and  |l0)  =  |l}®  |o) ,  with  the  result  that  no 

antisymmetric  currents  or  fluxes  are  generated  by  the  qubits  when  Molmer-Sorensen 
gates  are  employed.  A  coupled  SQUID  gradiometer  that  is  antisymmetric  with  respect  to 
the  flux  threading  the  two  cells  of  a  qubit  will  not  register  any  coupled  flux  until  a 
Hadamard  transform  is  carried  out  on  the  physical  qubits,  immediately  prior  to  a 
measurement.  The  two  basic  states  are  then  transformed  into  ±  combination  states  that 
have  antisymmetric  ‘figure  8’  circulating  currents  whose  flux  is  then  measurable  by  the 
SQUID  gradiometer.  The  coupling  of  each  qubit  to  an  LC  resonant  bus  is  made  by 
means  of  a  loop  that  encircles  both  cells  of  the  qubit  symmetrically.  This  prevents  the 
direct  mutual  inductance  coupling  of  the  bus  to  both  antisymmetric  microwave  gate 
pulses  and  the  measuring  SQUID  gradiometer. 

Pairing  the  qubits  into  logical  qubits  also  enables  Molmer-Sorensen  gates  to  be 
used,  since  the  phase  errors  that  would  ensue  if  QC  were  based  on  physical  qubits  are 
cancelled  out  in  each  pair  of  qubits.  The  requirement  that  the  LC  resonant  loop  be 
smaller  in  extent  than  one  tenth  of  a  substrate  wavelength  to  remain  within  the  quasistatic 
limit,  constrains  the  number  of  qubits  that  can  interact  on  a  single  LC  resonant  bus  to 
around  N  —  ~50  -  100.  To  enlarge  the  system,  two  resonant  loop  buses  are  coupled  by 
means  of  a  ‘transfer’  logical  qubit  that  couples  concurrently  to  both  buses.  Direct 
interaction  of  the  loops  is  prevented  by  coupling  them  antisymmetrically.  With  more  than 
one  transfer  qubit  per  bus,  the  buses  may  be  coupled  into  a  scalable  open  network.  The 
use  of  Molmer  -Sorensen  gates  insures  that  the  quantum  states  of  the  buses  and  transfer 
qubits  will  be  returned  to  their  pre-gate  status.  Since  the  logical  qubits  do  not  evolve  in 
time,  it  should  be  possible  for  the  system  of  coupled  buses  and  qubits  to  maintain 
quantum  coherence. 
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